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Abstract. In this paper, we solve explicitly the optimal stopping problem with random 
discounting and an additive functional as cost of observations for a regular linear diffu- 
sion. We also extend the results to the class of one-sided regular Feller processes. This 
generalizes the result of Beibel and Lerche [4, 5] and Paulsen [13]. Our approach relies 
on a combination of techniques borrowed from potential theory and stochastic calculus. 
We illustrate our results by detailing some new examples ranging from linear diffusions 
to Markov processes of the spectrally negative type. 



1. Introduction 

Consider a one-dimensional regular diffusion X = (Xt)t>o with state space E = (/, r), an 
interval of M, defined on a filtered probability space (J^()j>o,P). We denote by {Fx)xeE 
the family of probability measures associated to the process X such that Pa;(-'^o = x) = 1, 
and by E^; the associated expectation operator. Next, let be the family of all stopping 
times with respect to the filtration F (= {J^t)t>o)- In this paper, we are first concerned 
with the study of the following optimal stopping problem, for any x G E, 

(1) sup E, [e-^^g{XT) - Ct] 

where 5 is a non- negative continuous function on ^ = {At)t>o is a continuous additive 
functional of the form 

(2) At= f a{X,)ds, 

Jo 

with a a continuous function on E such that a{x) > for all x £ E, and, for any t > 0, 



(3) Ct= / c(X,)e-^^ds, 

Jo 

with c a non- negative continuous function on E. We also aim to solve this problem in 
the case X is a Feller process of the spectrally negative type, that is when it has only 
negative jumps. It is of common use to name g as the reward function, C = {Ct)t>o as the 
cost of observations and A as the random discount factor. We mention that in the case 
^ = and c = the problem (1) has been studied by Dynkin [11] and Shiryaev [33] in the 
general framework of Markov processes. Moreover, the case c = 0, that is without cost of 
observations, has been intensively studied in the literature for one-dimensional diffusions. 
In particular, Salminen [31] by means of the Martin boundary, suggested a solution to this 
problem in terms of the excessive majorant function. More recently, Beibel and Lerche 
[4, 5], relying on martingales arguments, solve this optimal stopping problem explicitly. We 
also mention that, by using standard fluctuation theory, Kyprianou and Pistorius [18] offer 
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solution to some optimal stopping problem arising in financial mathematics for appropriate 
diffusions. A related result on optimal stopping problems for one-dimensional diffusions 
with discounting has been presented by Dayanik and Karatzas [10]. They characterize 
excessive functions via generalized concavity and determined the value function as the 
smallest concave majorant of the reward function. In the former case, the value function 
is given as the solution of a free boundary value problem associated to a second order 
differential operator which is the infinitesimal generator of the one-dimensional diffusion 
X. The term free, which comes from the a priori, unknown region where the problem 
is investigated, forces one to set up an artificial boundary condition of Neumann type to 
get a well-posed problem. This is the so-called smooth fit principle. All these techniques 
are well explained in the book of Peskir and Shiryaev [26]. The literature regarding 
optimal stopping problems associated to diffusion with jumps is more sparse and focused 
essentially on the study of specific examples. In this vein, we mention the paper of Alili 
and Kyprianou [1] where the authors deal with the issue of pricing perpetual American put 
options in a market driven by Levy processes. We also indicate that Baurdoux [3] solved 
an optimal stopping problem associated to generalized Ornstein-Uhlenbeck processes of 
the spectrally negative type. 

In this paper, we propose to solve the optimal stopping problem (1) with a cost of 
observations of the form (3) for one-dimensional regular diffusions. Our strategy can be 
described as follows. First, by a time change device, we reduce the optimal stopping 
problem with random discounting to a one with a deterministic discount factor but asso- 
ciated to an appropriate time change diffusion. Then, by an argument of potential theory, 
we transform the problem (1) to an optimal stopping problem of the same form with a 
new reward function but without cost of observations. We proceed by using a result of 
Shiryaev [33] which states that in our context the optimal stopping time is the first exit 
time of the process from a Borel set. Finally, with this information at hand, we can use a 
Doob's /i-transform technique, with a proper choice of the excessive function, to transform 
our problem to an optimization problem which has been studied in details by Beibel and 
Lerche [5]. 

The remaining part of the paper is organized as follows. In the next section, we overview 
the basic facts about one-dimensional diffusions. In Section 3, we state and proof our main 
result which consists on solving the general optimal stopping problem (1). We also show, in 
Section 4, how to generalize our result to the class of regular one-sided Feller processes. The 
last section is devoted to the treatment of new examples ranging from linear diffusions to 
processes with one sided jumps such as spectrally negative Levy processes and self-similar 
positive Markov processes of the spectrally negative type. 



2. Preliminaries 

In this part, we provide some well-known facts about linear diffusions which can be 
found for instance in Ito and McKean [14] and in Borodin and Salminen [8]. We recall that 
(0, (J^t)t>o,P) is a filtered probability space. We consider a linear diffusion X = {Xt)t>o 
with state space E, as the solution to the stochastic differential equation (for short SDE) 



(4) 



dXt = b{Xt)dt + (j{Xt)dWt 
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where W = (Wt)t>o is a one-dimensional Brownian motion. It is supposed that a and b 
are continuous and cr{x) > for all x £ E. We assume that X is regular, that is 

p^.(rj, < +00) > 0, yx,y£E, 

where Ty = inf{i > 0;Xt = y}. The transition semigroup {Pt)t>o maps Ci,{E), the space 
of bounded and continuous functions on E, into itself. It follows that X is a Feller process. 
Moreover, for every t > and x € E the corresponding measure A 1 — > Pt{x, A), with A a 
Borel set, is absolutely continuous with respect to the speed measure m, a positive u-finite 
measure on E. More specifically, we have 

Pt{x,A)= / pt{x,y)m{dy) 
J A 

where Pt{-, ■) stands for the transition probability density which may be taken to be pos- 
itive, jointly continuous in all variables and symmetric. The scale function s of X is an 
increasing continuous function from E to M, satisfying, for any a < x < b, 

s{b) — s{a) 

and is given by 

s'{x) = exp (-2 / -^i^dz 



We also recall that the infinitesimal generator L of X is the second order differential 
operator given, for a function / G C^{E), the space of infinitely continuously differentiable 
functions with compact support, by 

Lf{x) = ^a\x)f"ix) + bix)nx). 

Next, from the general theory of one-dimensional diffusion (see [14, p. 128]), the Laplace 
transform of the first hitting time Ty is expressed, for any g > 0, as follows 



r h+{x) 

X <y, 



(5) E, [ 



h-{x) 



x > y. 



I hq (y) 

where h'^ (resp. h~) is the increasing (resp. decreasing) continuous solution to the differ- 
ential equation 

(6) Lh{x) = qh{x), x e E, 

with appropriate conditions at the non-singular boundary points. These functions h'^ and 
h~ are called the fundamental solutions of the equation (6). They are linearly independent 
and their Wronskian is defined by 

and the Wronskian with respect to d/ds{x) denoted by Wq is the constant given by 
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Moreover, for any q > 0, the Green function or the q-potential density u'^ is defined as the 
Laplace transform of the transition probabihty density, i.e. 

f + OO 



r+oo 

/ e-'^*pt{x,y)dt, Vx,?/ G E. 
Jo 



In particular, we have 



(8) u\x,y) 



Wg^h^{y)hg{x), x>y. 



We say that the process X is recurrent if and only if lim u'^(x,y) = 00, for all x,y £ E, 

which is equivalent to T'x{Ty < 00) = 1, for all x,y £ E. A diffusion which is not recurrent 
is called transient. In this case, the potential u defined as 

u(x,y) = lim u'^(x,y) 

is finite for all x,y £ E. Finally, we mention that if X is transient with lim X^ = r then 
for any x < y 

r+oo 

u{x,y) = / pt{x,y)dt, 
Jo 

= s{r) - s{y). 

3. Optimal stopping problem for linear diffusions 

Our aim is now to find the value of the function V^r- defined as the solution of the 
optimal stopping problem (1), that is 

V/,(x) = sup [e-^^ff(XT) - Ct] 

where A = {At)t>Q is a continuous additive functional of the form (2) and the cost of 
observations Ct = Jq c(Xs)e~"^*(is with c and g non-negative continuous functions. Our 
main result is stated in Theorem 3.9 below. It consists on reducing the optimal stopping 
problem (1) into a new one which can be described as follows. On the one hand, it has 
a modified reward function but without both cost of observations and discounting factor. 
On the other hand, it is associated to a new diffusion obtained from the original one by 
a random time change and by a Doob's /i-transform. It turns out that solving this latter 
optimal stopping problem amounts to finding the solution of an optimization problem 
which has been studied by Beibel and Lerche [5]. More precisely, our approach can then 
be split into the following three steps. 

1. First we time change X by the inverse of the continuous increasing functional A and 
we use the well-known fact that in our context the two processes have identical hitting 
time distributions. Hence, we may consider without loss of generality the problem (1) 
with linear discounting, i.e. At = qt for some constant q > 0. 

2. Then, we characterize the potential associated to the functional C and we show how to 
reduce our problem to an optimal stopping problem without cost of observations but 
involving a new reward function. 
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3. Finally, we borrow an idea of Williams [34] and Pitman and Yor [27] for constructing 
conditioned diffusions by the method of /i-transform. We transform the problem de- 
scribed in the item 2 to an optimal stopping problem without discounting factor which 
has been solved by Shiryaev [33]. 

3.1. Time change for nonnegative additive functional. We start our program by 
considering that ^ is a non-negative continuous additive functional of X of the form (2) 
and we assume, without loss of generality, that A^o = oo P^; a.s. Since A is a continuous 
increasing function, it admits an inverse functional which we denote by V and given, for 
ah t > 0, by 

Vt = inf{s >0;As>t}, 



-ds 



where Yt = Xvt for any t > 0. Moreover, if X is the solution to the SDE (4), then it 
follows, from the Ito's formula, that Y is the unique solution to the SDE 

dYt = (-] {Yt)dt + (^) iYt)dWt 



where is a Brownian motion with respect to the filtration = {J-Vt)t>o- The process 
Y remains a linear diffusion with respect to the filtration F^. In particular, it is a Feller 
process, see e.g. Lamperti [19], and its infinitesimal generator takes the form 

(9) L^/(x) = ^L/(x) 

for a smooth function / on E. Next, we consider an open interval B G E and denote 
by the first exit time of the process Y from B, it is plain that we have the following 
identity 

= Arpx a.s. 
We are now ready to state the following. 

Lemma 3.1. For any x £ E, we have, with the obvious notation, 

(10) sup [e-'^^giXr) - Ct] = sup E, U-^giYr) - C (-) {Ys)e-'ds 
where Y is characterized by its infinitesimal generator (9). 

Proof. From [30, III. 21], we have that for every F^-stopping time T, At is an F^- stopping 
time. Thus, we obtain that 



sup E, [e-^^g(XT) - Ct] = sup E, e-^g{Xvs)- [ ' c{X,)e 

SG^L I Jo 



dv 



The proof follows by performing the change of variable u = A^ in the integral on the 
right-hand side of the previous identity. □ 

Consequently, in the sequel, we can assume, without loss of generality, that the additive 
functional A is linear, i.e. At = qt for some q > and the cost of observations is Ct = 
J^c{Xs)e-^'ds. 
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3.2. Get rid of the cost of observations. Let us now introduce the function 5 defined, 
for any x £ E, hy 

(11) S{x) = [Coo] . 

In the following, we provide an expression of 5 in terms of the characteristics of X and 
give some conditions under which it is continuous and finite. 

Lemma 3.2. For any q > and x £ E, we have 

6{x) = I hgi^x) / {y)c{y)m{dy) + h't{x) / h' {y)c{y)m{dy) 



where Wq stands for the Wronskian of and with respect to the scale function s, as 
defined in (7) and m is the speed measure of X. Moreover, if c satisfies the integrability 
condition for any x £ E, 



(12) hq{xV y)hj{x A y)ciy)m{dy) <oo, 

then 6 is continuous and finite on E. 

Proof. By using Fubini's Theorem, we obtain that 

poo r 

5{x)= / e-i'E,[c{Xs)]ds= / u'i{x,y)ciy)midy) 



hq{xV y)h^{x Ay)c{y)m{dy) 



E 



where we have used the identity (8). The proof of the claims follows readily. □ 
Remark 3.3. We note that if g = and X is transient with lim Xt = r then 5 is given by 



t—>-oo 



(13) 6{x) = / {s{r) - s{y))c{y)m{dy). 

Je 

In this case, we mention that Khoshnevisan et al. [15] identify the law of the perpetual 
integral functional Coo of a transient diffusion as the law of the first hitting time of a 
random time change diffusion. 

We are now ready to state the following. 

Lemma 3.4. // 6 is finite on E then, for any x £ E, we have 

sup E, [e-^^giXr] - Ct] = sup E, [e"^^ {giXr) + 6{Xt))] - 6{x). 

Proof. Note that for any F-stopping time T, we have the identity in law 

Coo = Ct + e-'i^Coo o 0t 

where {6t)t>o stands for the shift operator, i.e. for any t,s >0, 6tw{s) = w{t + s). Since 
6 is finite, the strong Markov property yields 



sup E, [e-'^'giXT) - Ct] = sup E, [e'^^ {g{XT) + S{Xt))] - S{x). 



□ 
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A nice consequence of the previous result is that the general optimal stopping problem 
(1) is equivalent to an optimal stopping problem without cost of observations. Before 
stating our next result, let us introduce a few further notation. Let O denotes all the open 
subsets of E containing the starting point x of X. Let T,q be the class of stopping times 
of the form Tb = inf{t > 0; Xt B} where B £ O. 

Lemma 3.5. Let f be a continuous non-negative function. Then, for any q > and 

X £ E, 



(14) 

Proof. Let X be defined by 



sup [e-'i^f{XT)] = sup [e-''^^ f [XtJ] . 



Xt 



Xt 

d 



if 
if 



t < e, 



t > en 



where e^ is an exponential variable of parameter q > taken independent of F and d is 
a cemetery state. Note that X is always transient and clearly with a function / as above 
and using the convention f{d) = 0, we have, for any x £ E 



E, 



fiXt 



fiXt)] 



Therefore, there is a one-to-one correspondence between the excessive functions for X and 
the g-excessive ones for X (see Definition 3.6 below). Moreover, the g-excessive functions 
of the Feller process X are lower semi-continuous. Then, from Shiryaev [33, Cor. 3 p. 129], 
we obtain that 



sup E^. 



fix. 



T 



sup E^. 



f{XT,) 



Hence 



sup E, [e-^^/(XT) 



sup E, [e-«^s/(XT5)] 



□ 



3.3. Doob's h-transform. Our aim in this part is to show how to transform an opti- 
mal stopping problem with discounting factor to an optimal stopping problem without 
discounting. To this end, we recall some basic facts on excessive functions and Doob's 
/i-transform and we refer to the book of Borodin and Salminen [8, II. 5]. 

Definition 3.6. A non-negative measurable function h : E M U {oo} is called q- 
excessive, q > 0, for the process X if the following two statements hold true: for any 
X £ E, 

(i) e'^'E^ihiXt)] < h{x), t > 0, 
(ii) hme-'^*EJh(Xt)] = h(x). 

A q-excessive function is called q-invariant if for any x £ E and t > we have 

e-^'E^ihiXt)] = h{x). 
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A function h is (^-excessive (resp. (/-invariant) if and only if the process e^''^h{Xf) is 
a positive supermartingale (resp. martingale). For every y S (1,^) the functions x i— )• 
u'^{x,y), X I—)- hq{x) and x i— >• h~{x) are g-excessive. These functions are minimal in the 
sense that any other arbitrary non-trivial g-excessive function h can be expressed as a 
linear combination of them. 

Definition 3.7. Let q > and h be a q-excessive function. For any x & E such that 
< h{x) < +00, and t > 0, we define the new probability measure as follows 

dF'i = e-'i'^$^dF^ on Tf 
h[x) 

The process X under the probability measure is called the Doob's h- transform (or q- 
excessive transform) of X. It is also a regular diffusion process and thus a Feller process. 

Next, we use an idea of Williams [34] and Pitman and Yor [27] for constructing con- 
ditioned diffusions by the method of /i-transform by means of the Laplace transform of 
first passage times. In fact, we slightly generalize their methodology by considering as 
g-excessive function the Laplace transform of the first exit time of an open set by X. To 
this end, let B ^ O and we recall that we denote by the first exit time from B hy X, 
i.e. 

TB = inf{t >0; XtiB]. 
For any x £ E, we write the Laplace transform of the stopping time Tb as follows 

Thus, without loss of generality, the continuity of X allows us to restrict O to open 
intervals (a, b) for some I < a < b < r. It is well known that the function cj)^ is solution 
to the following Sturm-Liouville boundary value problem 

{Ln(x) = qu{x), x G (a, 6), 
u{a) = u{b) = 1. 

In other words, cf)^ can be written as a linear combination of the fundamental solutions 
and h'^ . Setting 

it is plain that the mapping is a (/-excessive function for X. Thus, as in the Defini- 
tion 3.7, we define the probability measure as follows 

(15) drf = e-'i^Bh^{XTs)dF^, on Tr^+. 

The diffusion X under the family of probability measures P'*^ = {¥ll:^)x£E is transient 
with 

(16) ¥fiX^-GdB) = l, 

except if g = 0, where ^ stands for the lifetime of X under P'^^. Clearly, the probability 
P^ (^ < cxd) is either 1 or for all x. Moreover, since X is solution to the SDE (4), then 
under the probability P'^ , the diffusion X can be characterized as the solution of the SDE 

dXt = {b{Xt) + log' {h^iXt)) a\Xt)) dt + a{Xt)dWt 
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where is a standard Brownian motion under the new probabihty P^^ 



Remark 3.8. As explained in [27], we take to be defined by the requirement that for 
each X E: B the process X runs up to the time Tb has the same law under P^ as it does 
under P^^ conditional on Tg < where is an independent exponentially distributed 
random variable with parameter q > 0. 

We are now ready to state and proof the main theorem of this section. 

Theorem 3.9. If 6 is finite then solving the problem (1) amounts to solving the following 
optimal stopping problem 

g{XTB) + KXTB)' 



(17) sup 



B 



h^iXTB 



where the probability P^^ is defined in (15). Moreover, there exists an open interval B* 
of E such that the value function of the optimal stopping problem (17) is 

^ ' hB'{u*) Beo/gas hB{u) 'I I > 

and the optimal stopping time is Tb* ■ 

Remark 3.10. If B* is a bounded interval, that is B* = {u\,ul), (18) reads 

g{u\) + 6{ul) ^ g{ul)+5{ul) g{u) + 5{u) 

hB'{u\) h^'iul) TeOu&dB h^iu) 

Proof. First, we deal with the additive functional C. Since 5 is finite, we have, from 
Lemma 3.4, that solving the problem (1) is equivalent to solving the following optimal 
stopping problem without cost of observations 

sup [e-"^ (giXr) + S{Xt))] . 

Then, from Lemma 3.5, we deduce that 

sup E, [e-^^ {g{XT) + S{Xt))] = sup E, [e'^^^ (^(XtJ + S{Xtb))] • 

Next, we use the Doob's /i-transform device. Let P^^ be the probability measure defined 
in (15), then we have 

-g{XTB)+6{XTB) 



(19) sup E, [e-'?^s(g(XTj + 5(Xr^))] = sup E^ 



h^{XTB 



which completes the proof of the first assertion. Finally, since X under P^ is transient, it 
is stated above that Tb < 00 a.s., the value function of the last optimal stopping problem 
is the solution to the following optimization problem 

g{u) +6{u) 



sup sup „ 

BdOu&dB fi [u) 



□ 



10 



CISSE, PATIE, AND TANRE 



To be complete, we mention that the optimization problem (19) has been studied in 
details by Beibel and Lerche [5]. We refer to their paper for more precise information 
concerning its solution for all possible choice of the reward function g. We also point 
out that, in the specific case a = and X is a standard Brownian motion, a similar 
optimization problem was studied by Graversen and Peskir [12]. 



4. Extension to one-sided regular Feller processes 

Let now consider X to be the cadlag modification of a one-sided regular Feller process 
defined on a filtered probability space (fi, (J^t)t>o, P) and taking values in an interval 
E G M. It means that X is a regular Feller process having jumps only in one direction 
which we assume, without loss of generality, to be of the spectrally negative type. That 
is X does not have positive jumps, i.e. Fx (sup4>o(Xt — Xt-) > O) = 0, Vx G E. For sake 
of simplicity, we also assume that the process X has infinite lifetime. We wish to extend 
the results of the previous section to this class of stochastic processes. In comparison to 
the diffusion case, the difficulty is that we do not have, in general, any information on the 
excessive functions for this class of Markov processes. Indeed, the infinitesimal generator 
associated to X is an integro-differential linear operator for which there does not exist 
general results regarding the solutions to the boundary value problem (6). Nevertheless, 
as explained in the following, the one-sided feature of X allows to identify an increasing 
excessive functions. 

Proposition 4.1. For any q > 0, there exists a unique increasing function hq'.E^ 
[0,oo], lower semi- continuous, such that, for any x,y £ E with x <y. 



[e 



-qTyl (^) 



In the case X is recurrent, the function x i— ^ h'^(x) is continuous. 

Proof. As a consequence of the regularity assumption, it is well known, see e.g. [7], that 
for each singleton {y} G X admits a local time at y, which we denote by = {V()t>o. 
The continuous additive functional is determined by its g-potential, liq, which is finite 
for any g > and given by 



oo 







From the definition of and the strong Markov property, we obtain the identity, see [7, 
Chap. V.3], 

uq[y) 

Next following for instance Ito and McKean [14] or Pitman and Yor [27], we write (pqix, y) = 
Ecc [e"''"^^] and for a fixed zq e E, we define 
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Next, using the fact that X has no positive jumps, we get for any x < z < y, and by 
means of the strong Markov property, 

(t>q{x,V) = (j)q{x,z)(j)q{z,y). 

Thus, from the identity 

4>q{x,z) = h^{x)/h^{z) 

we deduce that the choice of the reference point zq affects /i^ only by a constant factor. 
The monotonicity of the mapping follows readily from its definition and the absence of 
positive jumps for X. For the last assertion, we recall that x \-^E,x [e"''"^^] is a g-excessive 
function, see e.g. [7, p. 74]. By linearity, the mapping x i— )• h^{x) is also (/-excessive and 
thus finely continuous. Then, the claim is a straightforward consequence of the fact that 
if X is also recurrent then the fine topology coincides with the initial topology of E, see 
e.g. [2, p.243]. □ 

We point out that, in Patie and Vigon [25], by means of arguments based on the Martin 
boundary, Proposition 4.1 is extended to a larger class of homogenous Markov processes 
with only negative jumps. Following the proof of the Theorem 3.9, and observing that 
h^~°°''^\u) = {u) / {x) , the proof of the theorem below goes through verbatim. 

Theorem 4.2. With the notation used in Theorem 3.9 we assume that 

n* g{u) + 5{u) g{u) + 5{u) 

D := sup — — — — = sup — — — — < oo. 

u&E hg[U) u>x,u£E hg[U) 

Then, for any x G E 

sup E,. [e-i^g{XT) - Ct] = D*h+{x) - 6{x). 
If there exists a point u* > x such that D* = ^ , L , then the optimal stopping time 

hq{u ) 

of the problem (1) is given by 

=mi{t >0, Xt> u*}. 

Remark 4.3. Since Lemma 3.1 applies also in this more general framework, one could have 
also considered a random discounting factor in the previous result. 

5. Examples 

We now illustrate our methodology by presenting the solutions to some new optimal 
stopping problems. We consider both the diffusion case and also the case when the pro- 
cesses are of the spectrally negative type. 

5.1. An optimal stopping problem with cost of observations for one-sided Levy 
processes. Let Z = {Zt)t>Q be a spectrally negative Levy process starting from x G M, 
that is a process with stationary and independent increments having only negative jumps. 
Plainly, the law of Z is characterized by ■i/'j the Laplace exponent of Zi, which admits the 
following Levy-Khintchine representation, for any n > 0, 



(20) V'(^i) = ^^^^^ + bu+ I (e"^ - 1 - 'uxl{_i<2.<o})i^((i: 







IX 
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where 6 G M and o" > and the measure v is such that / (1 A y^)u{dy) < oo. Next, 

J —oo 

recalhng that ip is a convex function on [0, oo) with lim^^oo = +oo, we denote by 9 the 
non-negative largest root of the equation ip{u) = 0. We also mention that being continuous 
and increasing on [6,oo), ip has a well-defined inverse function : [0,oo) — t- [^,oo) which 
is also continuous and increasing. We refer to the excellent monographs of Bertoin [6] and 
Kyprianou [17] for background on Levy processes. We now consider a perpetual American 
option in a market driven by Z where the agent takes into account some costs from hedging, 
which might come from some transaction costs or liquidity issues. We assume that the 
agent hedges continuously and that he chooses the cost of observations c(x) = exp(7j;) 
to be of exponential form. The payoff at time t > of such a product can be written as 
follows 

(21) e-^'gie^')- [' e'^^e"''' ds 

Jo 

where g > is the risk-free rate, 7 > — 1 and 5 is a smooth function. Next, assuming that 
q = it is easy to check that one may choose the probability as the risk neutral 
probability measure. Therefore, in the sequel, we suppose that the characteristics b, a and 
are chosen such that 

(22) g = V'(i)- 

5.1.1. The Brownian motion with drift case. We start with the case where Zf = W^^^ = 
bt + Wt is a Brownian Motion with drift b starting from and the reward function g is 
defined, for any < L < K, hy 

' L — X, if X < L, 

0, if L < X < K, 
X- K, if x>K. 



(23) g{x) 



This function corresponds to the payoff function of a strangle option which is a combination 
of a put with exercise price L and a call with exercise price K. In this case, the condition 

(22) is fulfilled if g = — + 6. We want to compute the constant 

m Tr rl LI ] / 

Vw = sup Eq 



Jo 



The case without cost of observations, i.e. c = 0, has already been studied by Beibel and 
Lerche [4]. Next, it is well known that the functions /i^ and h~ defined in (5) are given 
by 

h+{x) = Die"'^ h-{x) = Dse"^^ 
where ai = —b+ ^J2q + b'^ and a2 = —b — \/2g + 5^, and Di, D2 are positive real numbers. 
The function 6, defined in (11), is finite if g > 67 + — and is given by 

<5(^) = 7 ^e^- 

q — 'jb — 7/2 
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Next, we introduce the function 

Gp{x) 



g{x) + 6{x) 



which, according to Theorem 3.9, give the solution of our problem. Then, if g > 76 + 



7 



and b > —1 then qi > 1 and 02 < —1. Thus, we verify easily the following inequalities: 
(24) sup (e-"i^(5(x) + g{x))) > sup (e-°i^(,5(x) + g{x))) > 



a;<0 



x>0 



and 
(25) 

Note that 



sup (e-°2^(,5(x) + g{x))) > sup (e-"2^(,5(x) + g{x))) > 0. 



x>0 



x<0 



and 



supGp(x) = sup Gp{x) 

x>0 x>0; 5{x)+g{x)>0 



supGp(x)= sup Gp{x) 

x<0 x<0; 5{x)+g{x)>0 



Ct2X 



inf 



x>0; S{x)+gix)>0 g{x) + 6{x) 



inf 



pe 



x<0; S{x)+g{x)>0 g{x) + 6{x) 



Then for all p G (0, 1) 



1 



< sup Gp{x) < - sup (e"°i^(g(x) + S{x))) < +00. 

x>0; S{x)+g{x)>0 P x>0 



and 



1 



< sup Gp{x) < sup (e-"2^(5(x) + 6{x))) < +00. 

x<0; 5{x)+g{x)>0 ^ — P x<0 

We assume that log(L) < = < log(i^), q > 1/2 + 6 and 6 > -1. Then (24) 
and (25) hold and as in Beibel and Lerche [4, Lemma 1, p. 98] , there exists a number 
p* e (0, 1) such that 

sup Gp* (x) = sup Gp* (x). 

x>0 x<0 

Let xi,X2 and p* be solutions with xi > log K, X2 < log L and p* G (0, 1) of the 
following system 

( {q-lb- 7V2)(e^i -K) + e-^i {q - lb - 7V2)(L - e^^) + 



(26) 



-X2 



Let M* 



p*f.aixi _p*ya2Xi 

6^^76 + 772 -g) 

(q-lb- 7^/2) (L - e^2) + g-^a 

6^1(^-76-7^/2) -e-^i 
I (e^i - K){q --fb- 7V2) + e"^! 

{q-lb- 7V2)(e^i - -fC) + e'^i 



p*e"i^'i + (1 -p*)e"2a:i 



p*aie"i^2 ^ (1 -p*)a2e"2'^2 

p*e"l^2 + (1 -p*ya2X2 

p* aie^i^i + (1 - p*)a2e"2'^i 
p*e°^i^i + (1 - y>*)e"2^i 
(g-76-772)(L-e^2) + eX2 
p*e°i^2 + (1 _p*)e"2a;2 



. Then, 



sup Eq 



5(^: 



M* 



76 - 72/2 ' 
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and, the optimal stopping time is 

T(,^,,,) =inf{t >0; W^f) i (xi,X2)}. 

5.1.2. The spectrally negative Levy case. In the second example, we consider as the reward 
function 

g{x) = {x-K)+. 

We write, for any x £ E, 



(27) Vz{x) = sup 



:)+-«/" 

Jo 







where a > and we recall that we choose q = V'(l)- 
Proposition 5.1. We assume that = V'(l) ~ V'(7) > ^.''^d 

(28) x<x* = -logi 



7 V(l-7)a 
T/ien 

V^(x) = e^-^ (e^ -i^)+ + a^ - a— 

V J 

and the optimal stopping time is Tx* = inf{t > 0; Zt > x*} . 

Proof. First, by means of Fubini's theorem and using the fact that p^ > 0, we get easily 
that, for any x £ E, 

roo poo 'yx 

5{x) = aEx / e-«*e^^*dt = ae^^ / e-^^+^^^^^dt = a—. 

Jo Jo P-y 

Then, from Lemma 3.4, we deduce that 

Vz{x) = sup [e-"^ ((e^^ - K)+ + 6{Zt))] - S{x). 

Next, we recall from Bertoin [6, p. 189-190] that, for any x < y, 

Ex [e"''^^] =e-<^M(f~^-), r>0. 

Then, writing 

^ I git— a; ' 

we have that 

G'(u) = e^-" (i^l{„>iogX} + ^^^^e^") ■ 
Since G'{u) > on (— oo,x*] and G'{u) < otherwise, we deduce from (28) that 

supG(n) = supG(n) < oo. 

u>x 

The proof of the claim follows then by applying Theorem 4.2. □ 
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5.2. Optimal stopping problems associated to self-similar positive Markov pro- 
cesses of the spectrally negative type. Let Z = (Zt)t>o be again a spectrally negative 
Levy process. We introduce, for any a > 0, the process X = {Xt)t>Q defined by 

(29) log{Xt) = ZA„t>0, 

where 



At = inf{s > 0; = [ e"^'' dr > t}. 
Jo 



We denote its law by Q^^ when it starts from x > 0. Lamperti [20] showed that X is 
an a-self-similar positive Markov process on (0, oo), i.e. a Feller process which enjoys the 
following a-self-similarity property, for any c > 0, 

(30) {{Xco^t)t>o,Qx) = {icXt)t>o,Q.) . 

It is plain that X is also of the spectrally negative type, in the sense that it has no positive 
jumps. Next, we recall that the law of Z is characterized by its Laplace exponent, ip, which 
is of the form (20). In the sequel, writing 9 for the largest root in [0, oo) of the equation 
1/^(9) = 0, we assume that the following conditions 

(31) 9 <a and lim = oo 

hold. The first condition secures that the lifetime of X is infinite since in the case < 9 < 
a, we consider X to be the unique recurrent extension which hits and leaves continuously, 
see Rivero [29] for more details. Under the second condition, the paths of the process X 
are of unbounded variation on any compact interval and the process X is regular. Next, 
we introduce more notation taken from Patie [24]. Define for any integers n 

n 

a„(V',a)""^ = V'(afc), ao(V',") = l, 

k=l 

and we introduce the entire function X which admits the series representation 



n=0 

It is important to note that whenever 9, the largest root of the equation ip{9) = 0, satisfies 
< a, it follows that all of the coefficients in the definition of T^^a are strictly positive. 
Then, Patie [24, Theorem 2.1] characterized the Laplace transform of 

Ta = inf{t > 0; > a} 

as follows. Suppose that < x < a. Then, for any q > 0, we have 

Next, we introduce the Ornstein-Uhlenbeck process associated to X which is defined, for 
any t > 0, by 

(33) Ut = e^{-t)X,^^t),t>Q. 
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where ex{t) 



X 



a\ and we write vx{t) 



log(l+At) 



the continuous increasing 



A ' ^ 

inverse function of e\. These processes were introduced and studied by Carmona et al. [9]. 

In particular, they proved that C/ is a Feller process and under the conditions (31), U has 

also infinite lifetime and is regular. In [22], the author computed the Laplace transform 

of the first passage times above of U as follows. With the obvious notation, for any r > 

and < .T < o. 



(34) 

where 

(35) 



—tl>,a 



2;/>,a('';a") 



E°° T{q + n) 
^^^7-^a„('0;a)x'" 



n=0 



and r stands for the Gamma function. By means of classical criteria on power series, it is 
easily seen that, under the second condition in (31), that I^^a{q;x) is an entire function 
in X and is analytic on the domain {q G C; ^{q) > —1}- We shall also need the following 
representation of the function I^^a{Q', x) 



(36) 



which is readily obtained by using the integral representation of the gamma function 
T{q) = e~'^r'^~^dr, ^{q) > 0. We postpone to the end of this section the description of 
some specific examples of these power series. Let us assume that g is a continuous function 
and q, /3 > 0. We are ready to introduce the following optimal stopping problems 



Vf(x)^ sup E^e-^^giXT) 



= sup E, [e-'^^g (Ut)] 



V,^'^(x) 



sup Ea; 



9{Ut) 



Vf '?(x) = sup E, 

TGEoo 



e-'^^'g 



. At-- 



f 

Jo 



U7"ds, 



l + /3/o e^^-ds^ 

We note that in the case Z is a Brownian motion with drift, the last optimal stopping prob- 
lem is intimately connected to the so-called integral option problem studied by Kramkov 
and Mordetski [16]. 

Proposition 5.2. Let us write for some function h 

a,(/i . = arg max --— • 

u>x,u€E n[u°') 

(1) If = a^^iX^^q.)) exists and x < a* then 

T^,a{qx°') 



Vf(x) 



-g{a*)- 
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(2) If = (lii,,a (x'-^')) ^^^•^^•^ '^'^'^ X < a=K then 



(3) If a. 



q . . 



q . 



exists and x < a* then 



g{a*). 



(4) If = (li(),a (x'-^')) ^^^'^^•^ '^''^^ X < a* then 



2-0, a ( ^ ) X^ 



/n a// the above cases the optimal stopping time is given by Ta^ . 

Proof. The first item follows readily from the identity (34) and Theorem 4.2. Next, let 



\ni{t > 0; = y(l + aXt)^]. 



The Mellin transform of the positive random variable has been computed by Patie, 
see [22, Theorem 2]. However, for sake of completeness, we provide a slightly different 
proof here which relies on a device introduced by Shepp [32]. In the proof of [24, Theorem 
1], it is shown that the mapping x i— ?• I^^Q,(rx") is a r-eigenfunction for the infinitesimal 
generator of X. Thus, by the Dynkin formula, using the fact that the function X^^^ is 
increasing and applying the dominated convergence theorem, we deduce that 



Integrating both sides of the previous identity by the measure e 
theorem and the change of variable u = r{T^^ + x); we get 



r X 



^dr, using Fubini's 



1^1 ~l~ X^j/,q) 



q . . 



ip,a \ ^ ) X^ 



where we have used (36). Moreover, from the definition of U (33), we observe the identity 



(37) 

Hence, we obtain that 



t: 



a.s. 



-qT^ 



■ip,a 



(f;xx-) 



We deduce the item (2) by an application of Theorem 4.2. We complete the proof of the 
proposition from [22, Corollary 3.2]. □ 
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Remark 5.3. We note that the identity (37) combined with the item (2) of the previous 
proposition allows to solve the following non homogeneous optimal stopping time problem 



Vf«(x) = sup 
Indeed, we easily deduce that 



_(1 + 2AT))^ 



where is characterized by 

aiu) 



= arg max 



u>x,ueE -T- , fs.-'x/iia] 



In what follows, we provide some examples of the power series and we refer to [24] for the 
description of additional examples. 

The modified Bessel functions. We consider Z to be a Brownian motion with drift v >Q, 
i.e. il^{u) = ^u^ + vu and we set a = 2. Its associated self-similar process is well known to 
be a Bessel process of index u. We have 

r(-i/ + i) 

Thus, we get 

^2Ax) = {x/2Y''^T{-v + 1)I_^ [V2i 

where lv{x) = Y1'^=q Ttir^^+n+i) stands for the modified Bessel function of index z/, see 
e.g. [21, Chap. 5] and 

/ 3.2 

where i^, x) = XlnLo (yrn-^" stands for the confluent hypergeometric function of the 
first kind, see e.g. [21, 9.9], and {q)n = ^Yin)* stands for the Pochhammer symbol. 

Some generalized Mittag-Lejfler functions. In [23], the author introduced a parametric 
family of one-sided Levy processes which are characterized by the following Laplace expo- 
nent, for any 1 < a < 2, and 7 > 1 — a, 

(38) = ((« + 7-l)a- (7 -!)«)• 

Its Levy measure is absolutely continuous with a density / given by 

g(o+7-l)s/ 

where C is a positive constant. We focus on the case 7 = 1 in (38). We have ipi{u) 
il^{u) = {u)a and 



1 r(a(n + l)) 

a.n[.W]OL) = — , ao = 1. 

r(a) 
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Thus, the power series can be written as fohows 

1-i},a{^;x) = r(a)A^^„(aa;) 

where Mc,p{x) = J2n=o r(an+/3) (^^^P- ^lA^) = T.n=o r(an+^) ) stands for the Mittag- 
Leffler function of parameters a, /? > (resp. of parameters a, I3,q > 0) introduced by 
Prabhakar [28]. 
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